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Геоінформаційна технологія формування 
кадастрів емісій та поглинань парникових газів 
у лісовому господарстві Тернопільської області 


Розроблено математичні моделі для просторового аналізу процесів поглинання вуглекислого газу 
в процесі росту дерев та процесів емісії парникових газів унаслідок вирубування лісів на теренах 
Тернопільської області. Створено геоінформаційну технологію, яка дає можливість формувати 
георозподілену базу даних, будувати просторові кадастри емісій і поглинань парникових газів та 
оцінювати вуглецевий баланс у лісовому господарстві Тернопільської області. 


Вступ 


Збереження стабільності кліматичної системи сьогодні розглядається як одна з 
найважливіших проблем, що постали перед людством в останні десятиріччя. Як відомо, 
зміна клімату відбувається через збільшення антропогенних викидів парникових газів, 
найбільшою мірою двоокису вуглецю (СО), в атмосферу (1-31. 

Лісові екосистеми відіграють роль потужних поглиначів та акумуляторів вуглецю 
з атмосфери. Запаси вуглекислого газу і темпи його депонування у стовбурній деревині 
залежать від продуктивності лісів, їхнього стану, породного складу та вікової структури, 
а також режиму господарювання 1, особливо, інтенсивності антропогенних та природних 
порушень (4). Виходячи з цього, метою цієї роботи є представлення математичного 
апарату та геоінформаційної технології для просторового аналізу вуглецевого балансу в 
лісовому господарстві Тернопільської області. 


Специфіка об'єкта і формування цифрової карти 


У теперішній час ліси займають близько 13,3 Уо території Тернопільської області (31. 
Через вирубування цінних порід дерев відбувається зміна видового складу лісів: зі 
зменшенням частки дуба 1 бука зростає частка граба, осики, берези тощо. На даний 
час найважливішими лісоутворюючими породами на Тернопільщині є: дуб - 50 У, 
граб - 13 Зо, сосна- 8,8 Ус, ялина - 6,5 У лісового фонду. Доповнюють його бук, 
ясен, береза, осика, тополя, вільха, клен, модрина та інші. 

Найбільші площі лісових масивів розташовані у західній частині області. Дещо 
менше - на півночі (особливо у Кременецьких горах). Значні лісові масиви є в долинах 
Серету, Збруча 1 Стрипи. Майже всі ділянки, які вкриті лісом, займають широ- 
колистяні ліси: грабові, грабово-дубові, букові та інші. Враховуючи ці особливості 
області, для здійснення аналізу вуглецевого балансу вибрано широколистяні та 
широколистяно-соснові типи лісу, які мають найбільший вплив на стік вуглецю. 
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Для аналізу вуглецевого балансу як вхідні дані використано цифрові карти 
потенційної та реальної рослинності регіону. Карта потенційної рослинності відображає, 
який тип лісу міг би рости на тій чи іншій території Тернопільської області за умов 
невтручання людини. Цифрова карта реального лісу відображає, на якій території в 
теперішній час росте ліс. На основі цих двох карт сформовано вихідну карту, яка 
відображає реальні території, зайняті в теперішній час лісом, з розбиттям на основні 
типи лісу (рис. 1). 
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Рисунок І - Карта лісів Тернопільської Рисунок 2 - Карта поглинутого вуглецю 
області (2009 рік, т/км:) 


Множина елементарних об'єктів для аналізу 
вуглецевого балансу 


Для просторової оцінки вуглецевого балансу лісів Тернопільської області 
необхідно насамперед сформувати елементарні об'єкти, в межах яких здійснюватиметься 
цей аналіз. За основу такого формування вибрано елементарні ділянки лісу, які 
обмежені границями адміністративних районів, а також межами ділянок з однаковим 
типом лісу, що росте на цій ділянці. Загальне число таких елементарних ділянок у 
межах області визначається формулою 


ме М, (1) 

де М - загальне число елементарних ділянок; г - порядковий номер адміністра- 
тивного району, г -1,...К; К - число районів в області; М, - число елементарних діля- 
нок в г-у районі. Таким чином, отримуємо множину елементарних об'єктів цифрової 


карти, які надалі використовуємо для оцінки вуглецевого балансу в межах кожного об'єкта. 
Загальне число таких отриманих елементарних об'єктів в області позначаємо як М. . 
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Інформаційна технологія 


Відповідно до розробленої інформаційної технології здійснюється ітераційний 
процес для всіх елементарних ділянок (п - 1,..., Мо ). На кожній ітерації виконуються 
такі операції: 

- аналіз елементарного об'єкта 1 визначення, яка частина якого типу лісу є на 
цій ділянці, оцінювання поглинутого вуглецю в процесі росту дерев на цій еле- 
ментарній ділянці (обчислення кількості вуглецю акумульованого у вирощеній 
стовбурній деревині); 

- оцінювання частки вилученого вуглецю з елементарної ділянки (обчислення 
кількості вуглецю у вирубаній деревині); 

- аналіз вуглецевого балансу в елементарному об'єкті (різниця акумульованого 
та вилученого вуглецю). 

Після завершення ітераційного процесу і отримання даних щодо вуглецевого 
балансу для всіх елементарних об'єктів цифрової карти здійснюється візуалізація 
результатів. Результати інвентаризації містяться в таблиці вихідних даних (у вигляді 
георозподіленої бази даних). По кожному стовпцю цієї таблиці можна побудувати 
окремий шар цифрової карти, тобто представити результат у вигляді окремої тематичної 
карти, що дає можливість краще сприймати і аналізувати отриману інформацію |6|, (71. 


Математична модель процесу нагромадження вуглецю 


Оцінювання поглинутого вуглецю т-м типом лісу в п-й елементарній ділянці 
здійснюється за формулою 
а срари - Кидній » (2) 
де К, - коефіцієнт, який відображає вміст вуглецю у деревині т -го типу лісу 
3 . чо . - . 
(Мг/м'), Р, - річний приріст стовбурної деревини для т-го типу лісу (м'/км:), 


- площа т-го типу лісу в п-й елементарній ділянці (км'). Коефіцієнт Ки Є 


пт 
різним для різних типів лісу і обчислюється за формулою 
Ки З пт? 3) 
де г, - Частка сухої речовини у вилученій деревині, с,, - частка вуглецю в сухій 
біомасі. 


Результати числового моделювання для Тернопільської області (рис. 2) показують, 
що найбільше вуглецю поглинається у Бережанському, Борщівському та Кременець- 
кому районах, а також значна частина поглинається лісами Гусятинського, Зборівського 
та Чортківського районів. Такі результати можна пояснити великою кількістю лісів на 
території цих районів, а також складом лісу, адже в основному в цих районах 
поширені хвойні ліси, які поглинають більше вуглецю, ніж широколистяні. 


Математична модель вилучення вуглецю 


Спершу визначаємо частку вирубаного лісу у кожному елементарному об'єкті. 
При цьому робиться припущення, що ліс вирубують рівномірно в межах адміні- 
стративного району і рівномірно по породах, тобто 


Р, за віро (4) 


де И, - відносна вирубка в г-му районі, а, - абсолютна вирубка в цьому 


й 
районі, 5;,, - площа, зайнята всіма породами лісу в г -му районі. Дані про абсолютну 
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вирубку лісу в кожному районі можна рівномірно «рознести» по всіх елементарних 
ділянках лісу в цьому районі: 
У, з 10004, 8, зт но (5) 


де /, - об'єм вилученої деревини з п-ї ділянки (м'), а, - дані про вирубку лісу 


в районі (тис. м'), 5, - площа п-ї елементарної ділянки в аналізованому районі (км'), 


Ж Й оз 2 
Зінт р - загальна площа лісу в г-му районі (км'), п, - число елементарних 


об'єктів лісу в межах /-го адміністративного району. Далі отримані дані використовуємо 
для обчислення вилученого вуглецю. 

Щоб обчислити кількість вилученого вуглецю (Мг) з кожної елементарної 
ділянки лісу, використовуємо наступну формулу: 


сити» (6) 


де К, - коефіцієнт вмісту вуглецю у вилученій деревині т-го типу, яка росте на 
цій елементарній ділянці (Мг/м'), У, - кількість вилученої деревини з п-ї ділянки (м'). 
Статистичні дані про промислове вирубування лісу по адміністративних районах 
можна отримати на основі статистичного довідника |8|. Для здійснення просторової 
інвентаризації використано такі параметри, як середня продуктивність лісу (або 
приріст), коефіцієнт, що визначає долю вуглецю в сухій біомасі, коефіцієнт, який 
відображає запас деревини на одиницю площі. Результати числового аналізу показу- 
ють, що, наприклад, у 2009 р. найбільше вилученого в результаті вирубування 
вуглецю було у Чортківському, Кременецькому 1 Бережанському районах (рис. 3). 
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Рисунок 3 - Вилучений вуглець (на рівні адміністративних районів, тис. т) 


Просторовий аналіз вуглецевого балансу 


На цьому етапі здійснюється оцінювання вуглецевого балансу для кожного 
елементарного об'єкта за формулою: 
А - Спр С, з» (7) 


де С, , - Кількість поглинутого вуглецю у п-му елементарному об'єкті, С, - 


пр 
кількість вилученого вуглецю. Вуглецевий баланс мп-ї ділянки оцінюємо тоді за 
допомогою формули 


по р о АННИ Я (8) 
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яка враховує приріст стовбурної деревини по всіх типах лісу на елементарній ділянці 
та вилучений у результаті вирубування вуглець. Підсумовування таких результатів 
по всіх елементарних ділянках дає вуглецевий баланс, відповідно, адміністративного 
району чи області в цілому. 

Розроблено геоінформаційну систему, яка відображає вуглецевий баланс на рівні 
проаналізованих елементарних об'єктів (як шар цифрової карти). На основі отриманих 
результатів просторової інвентаризації можна обчислити відповідні узагальнені пара- 
метри для Тернопільської області в цілому чи для адміністративних районів. Результати 
обчислень показують, що у Чортківському, Тернопільському та Кременецькому районах 
кількість вилученого в результаті вирубування вуглецю значно перевищує кількість по- 
глинутого в результаті росту дерев. Натомість, найкраща ситуація у Зборівському районі. 

На рис. 4 подано узагальнені дані для адміністративних районів. Наведені резуль- 
тати свідчать про те, що вилучений у результаті вирубування лісу вуглець у ряді районів 
Тернопільської області дещо перевищує вуглець, нагромаджений у процесі росту дерев. 
Саме тому вуглецевий баланс для цих адміністративних районів є від'ємним. Без- 
заперечними «лідерами» є Чортківський та Тернопільський райони області. Так, у 2009 
році вуглецевий баланс цих районів був від'ємним на Рани близько 12 тис. т. 
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Рисунок 4 - Карта нагромадженого та вилученого вуглецю на рівні 
адміністративних районів (2009 рік, тис. т) 


Аналіз динаміки вуглецевого балансу 
зо 


2000 2005 2006 2007 2008 2003 
Рисунок 5 - Сумарний вуглецевий баланс лісів Тернопільської області (верхня крива), 
сумарний вилучений вуглець в області (крива посередині) та баланс Чортківського 
району (нижня крива), тис. Т 
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Для аналізу динаміки вуглецевого балансу для кожного елементарного об'єкта 
і для адміністративних районів здійснено просторову інвентаризацію як акумульованого, 
так і вилученого вуглецю за ряд років, для яких були доступні статистичні дані. Резуль- 
тати обчислень подано на рис. 5 та рис. 6. Найкраща ситуація спостерігалася у 2000 р., 
а найгірша - у 2009 р. На рівні області найбільше вилученого вуглецю було у 2007 році. 


Висновки 


Опрацьовано підхід до формування елементарних об'єктів цифрової карти 
Тернопільської області, які відповідають окремим ділянкам та типам лісу і максимально 
пристосовані до реалізації просторового аналізу вуглецевого балансу в секторі лісового 
господарства області. 
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Рисунок 6 - Карта вуглецевого балансу на рівні адміністративних районів, тис. т 
(2000, 2005 -- 2009 роки, відповідно) 


Розроблена математична модель та алгоритм аналізу нагромадженого вуглецю 
дають можливість здійснювати просторову інвентаризацію акумульованого у стовбур- 
ній деревині вуглецю у процесі росту лісу (з урахуванням нових насаджень). Запропо- 
нована математична модель та алгоритм аналізу вилученого вуглецю дають змогу 
знаходити обсяг вуглецю, який втрачається в лісовому «резервуарі» у процесі проми- 
слового вирубування лісу. Розроблені математичні моделі відображають просторовий 
розподіл аналізованих величин і дають можливість знаходити вуглецевий баланс для 
елементарних об'єктів цифрової карти. 

На основі відомих статистичних даних здійснено числові експерименти, які 
показали, що процеси поглинання та вилучення вуглецю відбуваються нерівномірно 
на території Тернопільської області. Проілюстровані результати дають змогу побачити 
узагальнену картину вуглецевого балансу в області, виявити райони, у яких найгірша 
ситуація. Це дає можливість відповідним владним структурам приймати ефективні 
рішення щодо раціонального ведення лісового господарства. 
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Р.А. Бунь, Й.В. Слука 

Геоинформационная технология формирования кадастров змиссий и поглощений парниковьіх газов 
в лесном хозяйстве Тернопольской области 

Разработаньт математические модели для пространственного анализа процессов поглощения углекислого 
газа в процессе роста деревьев и процессов змиссий парниковьгх газов в результате вьтрубки лесов на 
территории Тернопольской области. Создана геоинформационная технология, позволяющая формировать 
геораспределенную базу данньх, строить пространственньте кадастрьт змиссий и поглощений парниковьгх 
газов, а также оценивать углеродньій баланс в лесном хозяйстве Тернопольской области. 
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